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Using 383 million BB pairs from the BABAR data sample, we report results for branching fractions 
of six charged B-meson decay modes, where a charged kaon recoils against a charmless resonance 
decaying to KK* or iq-K-R final states with mass in the range (1.2 — 1.8) GeV/c 2 . We observe a signif- 
icant enhancement at the low KK* invariant mass which is interpreted as B + — > 7/(1475)ii' + , find 
evidence for the decay B + — > rj(1295)K + , and place upper limits on the decays B + — > r](1405)K + , 
B+ -» /i(1285)if+, B+ -> /i(1420)Jf+, and B+ -» 0(1680)^"+. 

PACS numbers: 13.25.Hw, 14.40.Cs, 12.39.St, 12.39.Mk 



Charmless hadronic _B-meson decays have been of par- 
ticular interest due to their sensitivity to weak inter- 
action dynamics. The first observed gluonic-penguin - 
dominated decays, such as B — ► rf K and B — > 7r7^T [lj], 
allowed the study of CP violation in these decays with 
potential sensitivity to new physics @, Q . The relatively 
large B — ► 7/'ii" decay rate was also a topic of debate. 
However, little is known about the B meson decays to 
excited states of the 7/ and rj mesons. There are three 
candidates for the first excited states 77(1295), 7/(1405), 
and 7/(1475) and there is a possibility that they might 
include a gluonium admixture [B[ . This part of the pseu- 
doscalar meson spectrum remains uncertain after a few 
decades of studies [|1 0, H, 0, EH 11 1. A search for 
B-meson decays to these pseudoscalar states is the focus 
of this Letter. 

The r\ and rf candidates and their excited counter- 
parts, which we call generically rjx in this paper, have 
the quantum numbers J p = 0~ and decay strongly to 
at least three pseudoscalar mesons. Thus we look for 
the rjx — > KKtt and 777T7T final states. In the former 
case, the resonant structure KK* + KK* is of particu- 
lar interest and we refer to it as KK*. Previously, the 
K*K + K~ final state has been studied by BABAR inclu- 
sively 0. The J p = 1+ mesons A (1285) and /i(1420) 
and J p = 1" meson </>(1680) also appear in the mass 
range (1.2 — 1.8)GeV/c 2 in these final states. These 
resonances are considered in our search for the decays 
B + — > rjxK + and referred to by the generic nomen- 
clature r/x as well. Hermitian conjugation is implied 
throughout this paper unless stated otherwise. 

The B — > i]xK decay mechanism is expected to be 
dominated by the b — > s gluonic-loop penguin diagram, 
similar to the B — ► rj'K decay. The expected branch- 
ing fractions differ significantly depending on the rjx 
state [4j , following a pattern that early naive factorization 
models were unable to predict [lj|. The first attempt at 
unraveling the pattern in the branching fractions of B- 
meson decays with 77 and r\ 14J suggested including the 
interference within the quark flavor octet among other 
possible scenarios, but the predictions did not match the 
experimental data. More recent calculations find a larger 
predicted rate for B — > rj'K , in agreement with data, with 
inclusion of higher-order corrections [l5[ or "charming- 
penguin" contributions [l6j |: large theoretical uncertain- 
ties persist, partly due to insufficient experimental data. 
An admixture of a bound two-gluon state, gluonium, in 



r/x could also explain the enhancement of the branching 
fractions. 

Although the 7/(1295), 7/(1405), and 7/(1475) states 
are considered well-established [4(, their nature is still 
unknown. Partial wave analyses of the KKtt and 
777r7r spectra from past experiments, such as studies in 
Refs. 0, 1, ©, [To| . conclude that the meson spectrum in 
the (1-2— 1.8) GeV/c 2 range is described by a linear com- 
bination of the resonant states and a nonresonant phase- 
space contribution. The analyses in Refs. @, H| found 
that mass spectrum description without interference be- 
tween the resonant and nonresonant contributions is pre- 
ferred. Therefore, in our analysis we adopt the model 
of three spin-zero resonances 7/(1295), 7/(1405), 7/(1475), 
three spin-one resonances /i(1285), /i(1420), 0(1680), 
and a phase-space nonresonant contribution without in- 
terference with the above states. Only four resonances 
are considered in each final state, KK* or 7/7r7r, accord- 
ing to their dominant decay modes as discussed below. 

We use a sample of (383 ± 4) million T{AS) -> BB 
events collected with the BABAR detector [ljj at the PEP- 
II e + e~ asymmetric-energy storage rings with the e + e~ 
center-of-mass energy y/s = 10.58 GeV. Momenta of 
charged particles are measured in a tracking system con- 
sisting of a silicon vertex tracker with five double-sided 
layers and a 40-layer drift chamber, both within the 1.5- 
T magnetic field of a solenoid. Identification of charged 
particles is provided by measurements of the energy loss 
in the tracking devices and by a ring-imaging Chcrenkov 
detector. Photons are detected by a CsI(Tl) electromag- 
netic calorimeter. 

We search for B + — > r/xK + where rjx decays to KK* 
and7/7r + 7r-. We reconstruct ifif* -> K >0 K ± tt^, ( K >a -> 
Kg — > 7r + 7r~, and r\ — > 77. Isospin symmetry implies 
that the final states K°K~it + + K°K + it~ and 7/7r + 7r~ 
constitute two thirds of r/x — * KK* and r/x — > t]tttt, 
respectively. 

We identify B meson candidates using two kinematic 
variables: tties = \/ s /4 — p% and AE — \fs/2 — Eg, 
where (Eb,Pb) is the four-momentum of the B candi- 
date in the e + e~ center of mass frame. We require tties > 
5.25GeV/c 2 and \AE\ < 0.1 GeV. The requirements on 
the invariant masses are 1.35 < rn K ^, < 1.8 GeV/ c 2 , 
1.2 < m I)7rw < 1.5GeV/c 2 , \m^ — m K o\ < 12MeV/c 2 , 
and 510 < 77i 77 < 570MeV/c 2 . The rjx invariant mass 
range is chosen to include the broad spectrum of states 
without extending it above the charm background pro- 
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duction threshold. 

We require the photon energies be at least 100 MeV. 
For the Kg candidates, we require the cosine of the angle 
between the flight direction from the interaction point 
and the momentum direction to be greater than 0.995, 
and the measured proper decay time to be greater than 5 
times its uncertainty. In the r\x — > KK* +KK* — > KKir 
decay channel, we require the Kn or Kn invariant mass 
to satisfy 0.85 < m K7T < 0.95 GeV/c 2 for either A" ± 7r T or 

K °7rT combinations. 

We use the angle Qt between the -B-candidate thrust 
axis and that of the rest of the event, and a Fisher dis- 
criminant J~l to reject the dominant e + e~ — > quark- 
antiquark background [l8[ . Both variables are calculated 
in the e + e~ center-of-mass frame. The discriminant com- 
bines the polar angles of the £>-candidate momentum vec- 
tor and its thrust axis with respect to the beam axis, and 
two moments of the energy flow around the B-candidate 
thrust axis [lij ]. 

We suppress the background from i?-decays into states 
with D or cc mesons by applying vetos on the invariant 
masses of their decay products. The remaining back- 
ground (less than 10%) comes from random combina- 
tions of tracks from B decays, and from B + — > KK*K + . 
When more than one candidate is reconstructed, we se- 
lect the one with the lowest combined \ 2 01 the charged- 
track vertex fit and of the invariant mass of the K ° or 77 
candidate relative to the PDG values [4(. 

We define the helicity angle 9-u as the angle between 
the direction of the B meson and the normal vector to 
the rjx three-body decay plane in the r\x re st frame. The 
ideal distribution is uniform, 7i 2 , or (1 — H 2 ) for r\x with 
J p = 0~, 1~, or 1 + , respectively, where H — cos9 H . 
The observed angular distribution can be parameterized 
as a product of the ideal angular distribution for a given 
spin and parity multiplied by an empirical acceptance 
function parameterized as a polynomial P(|7i|). 

We use an unbinned, extended maximum-likelihood fit 
to extract the event yields rij and the parameters £ of 
the probability density functions (PDF) Vj. The index 
j represents six event categories used in our data model: 
the B + — > i]xK + signal (four categories in each of the 
two rjx decay channels as shown in Table [TJ) , combinato- 
rial background (mostly e + e~ — > qq production with a 
few percent admixture of misreconstructed -B-meson de- 
cays), and a possible background from B — > KK* K (in 
the i]x — > KK* channel) or other B backgrounds (in the 
Vx ??7T7r channel) . The likelihood Ci for each candidate 
i is defined as Ci — ^2jiijVj(xi,^), where the PDF is 
formed from the observables x = {toes, A_E, Tl-, W, to}. 
Here to is the invariant mass of the r\x candidate. 

We use a relativistic spin-J Breit-Wigner amplitude 
parameterization for the invariant mass of an rjx res- 
onance with the nominal mass and width parameters 
quoted in Table |U We model the decay kinematics as 
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FIG. 1: Projections for B + KK*K + (left column) and 
B + — > ri-KirK + (right column) of (a,b) tties, ( c ,d) AE, (e,f) m 
with a requirement applied on the signal-to-background prob- 
ability ratio calculated with all variables except the one be- 
ing plotted. The extended mass region in (f) includes the 
77' resonance as a crosscheck. The nominal region is shown 
in the inset. The solid (dashed) lines show the signal-plus- 
background (background) PDF projections. The dotted line 
shows the total PDF projection excluding the rj(1475)K + 
(left) or T7(1295)ii" + (right) final states. The dash-dotted lines 
indicate the nonresonant component. The long-dashed line in 
(e) represents the cross-check with the 77(1475) resonance mass 
(mo) and width (F) parameters unconstrained, both resulting 
in larger values. 

r\x — * KK* — > K Ktt and r\x — * ao(980)7r — > t/tttt. For 
the rjx — > KK* mode, the rjx invariant mass parameteri- 
zation is corrected for phase space of the B + — ► KK*K + 
decay and averaged over the K* — > Ktt invariant mass 
values. We ignore the interference between the overlap- 
ping resonances because it averages to zero for resonances 
with different quantum numbers or because these reso- 
nances have different final states, such as 77(1405) and 
?7(1475). The former decays mainly to ao(980)7r (or di- 
rect KKir) and the latter mainly to KK* (4J]. We also 
ignore the interference between the resonant and nonres- 
onant decays based on indications from previous stud- 
ies of rjx decays @, H| and due to potentially different 
three-body structure. This interference effect would only 
increase the significance estimate because the hypothesis 
of zero yield is not affected and the likelihood of the nom- 
inal fit could only improve. The significance is defined as 
the square root of the change in 2 In C when the yield is 
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TABLE I: Summary of results for the B — > r\x K+ process studied with six B-decay modes and eight decay channels with the 
signal resonance and nonresonant model discussed in text, where r\x — * KK* — > KsK ± n^ 1 in the upper part and rjx — * r/n + n~ 
in the lower part. The mass mo and width T of six rjx states are quoted Q] with errors in parentheses. The number of 
signal events n sig with significance of the observed signal in parentheses, the product of the branching fractions B and the 
corresponding daughter branching fractions, the B + — > /i(1285)_R" + branching fraction, the corresponding 90% C.L. upper 
limits, and selection efficiencies e obtained from MC simulation are shown. The systematic uncertainties are quoted last. 



77x KK* 
m /r [4], MeV 

90% C.L. 
B{B + -» rjxK + ) B{r,x -> KK* 
90% C.L. 

e(%) 



r7(1475) 
1476(4) /87(9) 
155l?i+^ (7.5a) 
192 

(13.8 



1-21+U 
-19 

< 

+1.8+1.0 \ 
-0.6) 

X 10 



1.7 

< 17 
8.8 ±0.1 



10" 

■6 



0(1680) 
1680(20)/150(50) 
17±| ± 7 
< 39 

< 3.4 x 10" 6 
9.0 ±0.2 



7^(1405) 
1409.8(2.5)/51. 1(3.4) 
-12+* ± 1 
< 12 

(-1.2i°;^±0.1)10- 6 
< 1.2 x 10 -6 
8.4 ±0.3 



/i(1420) 
1426.3(0.9) /54.9(2.6) 
36±i4 ± 7 
< 56 

(2.7±S;g±0.5) 10~ 6 
< 4.1 x 10" 6 
10.7 ±0.3 



m /r [4], MeV 

90% C.L. 
B(B+ -> rjxK + ) B{r) X -» t/tttt) 

90% C.L. 
B(5 -> /i(1285)if + ) 

90% C.L. 
e(%) 



77(1295) 
1294(4) /55(5) 
131±g| ± 10 (3.5a) 

< 179 
(2.9tg;? ±0.2) 10~ ( 
< 4.0 x 10~ 6 



17.6 ±0.3 



A (1285) 
1281.8(0.6)/24.2(1.1) 
-30±i9 ± 14 
< 30 

(-0.8±g:| ± 0.4) 10~ ( 

< 0.8 x 10" 6 
(-1.5+^ ± L2 ) 10-6 

< 2.0 x 10" 6 
14.1 ±0.9 



77(1405) 



< 54 

-o.3±g;| ± 0.1) 10" 

< 1.3 x 10" 6 



16.5 ± 1.2 



A (1420) 



49Z 



± 11 



< 99 

(1.4 ±1.0 ±0.3) 10~ 6 
< 2.9 x 10~ 6 



13.5 ±0.6 



constrained to zero in the likelihood C 

The signal PDF for a given candidate i is the prod- 
uct of the PDFs for each of the discriminating variables. 
The combinatorial background PDF is the product of the 
PDFs for independent variables. The signal and back- 
ground PDFs are illustrated in Fig. [T] We use a sum of 
Gaussian functions for the parameterization of the sig- 
nal PDFs for AE, Toes, and Tl- For the combinatorial 
background, we use polynomials, except for toes and Tl 
distributions, which are parameterized by an empirical 
phase-space function and by Gaussian functions, respec- 
tively. The nonresonant B — > KK*K background is pa- 
rameterized the same as signal, except for the quantity 
to, which is described by a phase-space function. 

The PDF parameters (£) of the combinatorial back- 
ground are left free to vary in the fit, except for the pa- 
rameters that describe Tl and the toes endpoint, which 
are fixed to the values extracted from the data sideband 
region (to ES < 5.27GeV/c 2 or \AE\ > 0.07 GeV). The 
PDF parameters for other event categories are taken from 
Monte Carlo (MC) simulation [l^| and adjusted with 
B — > Dir calibration data samples. We allow the yields to 
become negative as long as the total likelihood function 
remains positive in the allowed ranges of the observables. 
We study the goodness-of-fit and validate the fit proce- 
dure using MC simulation and generated samples. 

In Table U we present the results of the fit. We observe 
a large charmless contribution in the B + — > (KK*)K + 



decay with a significant enhancement at the low KK* 
invariant mass, which is interpreted as ?7(1475) — > KK* 
from the decay B + — ► r/(1475)K + . We also see ev- 
idence for a nonzero B + — > rj(l2%5)K + yield in the 
77(1295) — + 777T7T channel. The significances are more 
than 7.5 and 3.5 standard deviations, respectively, in- 
cluding systematic uncertainties. The significance of the 
B + — > 77(1295)X + yield is obtained in the fit when all 
yields are restricted to be positive, thus reducing the 
significance from the nominal fit. The significance is 
calculated within the model of resonant and nonreso- 
nant signal contributions discussed above and in earlier 
work [J, |7|, [8(. We quote 90% confidence level (C.L.) up- 
per limits, taken to be the values below which lies 90% of 
the total of the likelihood integral in the positive branch- 
ing fraction or yield region. 

We repeat the fit by varying the fixed parameters in 
£ within their uncertainties to obtain the associated sys- 
tematic uncertainties. The biases from the presence of 
fake combinations or other imperfections in the signal 
PDF model are estimated with MC simulation. Addi- 
tional systematic uncertainties originate from other po- 
tential B backgrounds, which we estimate can contribute 
at most a few events to the signal component. As a 
cross-check, we repeat the fit with the particle identifi- 
cation on the recoil kaon reversed in order to enhance 
the B + — ► rjx^ + topology by more than a factor of ten 
compared to the nominal reconstruction, and find no ev- 
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idence for such a decay. The systematic uncertainties 
in selection efficiencies are dominated by those in parti- 
cle identification, track finding, and K° and r] selection. 
Other systematic effects arise from event-selection crite- 
ria, and the estimation of the number of B mesons. 

The states 77(1475)^(1680), and /i(1420) are expected 
to decay into the KKtt final state through KK* |4j . We 
cross-check the KK* dominance by removing the Kit 
mass requirement and find consistent results. With the 
present dataset we are unable to resolve intermediate 
states in the r/mr modes, such as p (770) and 0^(980) 
resonances. 

In the projection plots in Fig. [TJ for illustration pur- 
poses, the signal fraction is enhanced with a requirement 
on the signal-to-background probability ratio, calculated 
with the plotted variable excluded. The m projection 
plot in Fig. [1] (e) implies a possible difference of the sig- 
nal resonance parameters from the assumed values. We 
repeat the fit with the 7/(1475) resonance parameters mo 
and r unconstrained while constraining other fit param- 
eters to the values from the nominal fit. We find the 
mo and T central values to be larger, but still consistent 
with the nominal values within statistical uncertainties 
(1482 ± lOMeV and 108 ± 20MeV, respectively). We 
also repeat the fit with the m range extended up to 2.5 
GeV/c 2 and find good extrapolation of the fit results in 
the full range, apart from the narrow charm production 
contribution just above the 1.8 GeV/c 2 threshold. 

In summary, we have measured product branching 
fractions B(B + — > rjxK + ) x B(r)x — > KK*,r)nn) for six 
-B-decay modes that have not been studied previously, 
where rjx stands for either 7/(1295), 7/(1405), 7/(1475), 
/i(1285), /i(1420), or 0(1680). We observe a significant 
enhancement at the low KK* invariant mass which is in- 
terpreted as B + — > ri(lA7b)K + and find evidence for the 
decay B + — > ?/(1295)if + . These decays could be used 
to either test weak dynamics in the predominant b — > s 
gluonic-loop penguin transition or study the rjx compo- 
sition, including potential gluonium admixture. 
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